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ABSTRACT 

An orbital solution for the MACHO*05:34:41. 3-69:31:39 eclipsing binary system is 
presented, based on the published light curve and spectral data obtained with the 
2.15-m telescope at CASLEO. Based on these spectroscopic observations, this system's 
binary components were classified as 03 IP and 06:V respectively. The radial velocity 
data along with the published light curve were analysed with the Wilson-Devinney 
code to derive the following masses and radii for the components of this system: 
Ml = 41 ± 1.2 M0, i?i = 9.6 ± 0.02 R©, M2 = 27 ± 1.2 Mq and R2 = 8.0 ± 0.05 Rq. 
The solution shows that the system is in overcontact as one would expect from the 
derived masses and the very short orbital period (^ 1.4 days). 
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1 INTRODUCTION 

One of the consequences of the huge amount of data pro- 
duced by the MACHO project was the discovery of more 
than six hundred new eclipsing binary stars in the Large 
Magellanic Cloud (Alcock et al. 1997). Given that the 
brightest among these systems lay within the scope of 
moderate size telescopes for intermediate dispersion spec- 
troscopy, we began to acquire spectra of those objects us- 
ing the 2.15 — m telescope at CASLEOn (San Juan, Ar- 
gentina), with the aim of producing determinations of phys- 
ical parameters for them. One of the observed stars, the 
so-called MACHO*05:34:41.3 -69:31:39 , located in the Hll 
region DEM 242 (Davies, Elliott & Meaburn 1976) turned 
out to be a double hned 03If*-|-06:V system, a spectral 
type we would not expect for such a short period binary. 
Apart from the general surveys as the Digitized Sky Survey, 
there are no previous observations of this star. 

There are an intrinsically small number of early O-type 
stars, and an even smaller number of them in eclipsing bi- 
nary systems. No 03 star is known to be an eclipsing binary 
to date. Besides, the masses predicted from comparison of 
the luminosities of these very hot stars with numerical mod- 
els are of the order of 100 Mq , a value much larger than 
that inferred from observations of binary stars of similar 
spectral types. For these reasons, each newly discovered 03 
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star presents an opportunity to provide new data that con- 
tributes to resolving this problem. 

Of course, it must be considered that numerical models 
deal with single and isolated stars, while the binary studied 
here is an extremely close system, with a period of ~ 1.4 days 
and undoubtedly subject to very strong interaction effects. 
Consequently, the masses derived here must be considered 
with caution when compared with evolutionary models. 



2 OBSERVATIONS AND DATA ANALYSIS 

2.1 Photometric and spectroscopic data 

The photometric light curve for MACHO*05:34:41.3- 
69:31:39 was published by Alcock et al. (1997,1997a). Data 
are available for two bands, namely V and R. The V mag- 
nitude at maximum is 13.54, while primary and secondary 
minima have V — 13.91 and 13.89 respectively. The V — R 
color of the system is —0.15 (in the Kron-Cousins photomet- 
ric system), with no noticeable changes with phase. Uncer- 
tainties given by Alcock et al. rise to 0.07 in V and 0.03 in 
V -R 

During an observing run in January 2000 with the 2.15— 
m telescope at CASLEO (San Juan, Argentina), N.I. Morrell 
acquired 8 spectra near the times of quadratures. 

A REOSC spectrograph was used in its single disper- 
sion mode, with a 600 grooves mm~^ grating and a Tek 
1024 X 1024 CCD as detector, this instrumental configura- 
tion giving a reciprocal dispersion of 1.8 A plxel"^ in the 
blue-yellow region of the spectrum. The spectrograms cov- 
ered the region from 3800 A to 5440 A, with a typical reso- 
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Figure 1. CCD spectrograms of MACHO*05:34:41. 3 -69:31:39 
obtained at CASLEO, corresponding to botli quadratures. 



lution (as measured from the He-Ar comparison lamp lines) 
of 2 pixels and the signal to noise ratio ranged from 50 to 
100. 

The usual sets of bias, dark and flat-field frames were 
acquired during each observing night. The spectroscopic ob- 
servations were processed and analysed with IRAF Pj rou- 
tines. In order to check the system's stability for radial ve- 
locity measurements, radial velocities for the nebular [Olll] 
emission lines in our spectra were derived, obtaining an aver- 
age of -1-258 kms"'^ with a standard deviation of 22 kms~^, 
a radial velocity that is well within the expected range for 
LMC objects. 

Just a quick look to the obtained spectrograms revealed 
that the primary star is an 03 If* presenting strong Nv 4603 
- 4620 A absorption lines and emissions of Niv, He ll and 
very weak Nlll. The spectrum of the secondary is hard to 
classify due to relatively low dispersion data, but considering 
those observations where the lines are better resolved, it was 
classified as 06:V. Both components were classifled follow- 
ing the criteria described in Walborn & Fitzpatrick (1990). 
Spectra of MACHO*05:34:41.3 -69:31:39 obtained at both 
quadratures are shown in Fig. 1. 

To obtain the radial velocities, the N v A4603 - 4620 
A absorptions and the Niv 4058 A emission were consid- 
ered as representative of the motion of the primary 03 If* 
component, while for the radial velocities of the secondary 
component several Hel absorptions were averaged. The mea- 
surements were performed through single Gaussian fitting 
with IRAF routines. 



2.2 Light and radial velocity curve solution 

The radial velocities derived for the quadratures (Table 1) 
together with the published photometric light curve were 
analysed by means of the Wilson-Devinney code (Wilson & 
Devinney 1971, Wilson 1990), that is very suitable to model 
binaries with Roche-lobe flUing stars. 

From Schmidt-Kaler (1982), a temperature of 5QQQ0 K 
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Table 1. Observed radial velocities of the eclipsing binary 
MACHO*05:34:41. 3 -69:31:39 



HJD 


2.5 X 10'' + 


1555.604 


1557.588 


1562.745 


1563.757 


1564.613 


1566.712 


1567.622 


1571.600 


800 


600 


'» 400 


E 200 


>" ° 


— ZOO 


-400 



phase 



t 



primary secondary 



0.799 521 -127 

0.211 -7 659 

0.882 447 -28 

0.603 428 — 

0.212 -47 — 

0.706 495 -165 

0.354 55 633 

0.185 -24 613 
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Figure 2. Observed and modelled radial velocity curve for 
MACHO*05:34:41. 3 -69:31:39 . Open circles correspond to the 
primary component and filled ones stand for the secondary. 



was assigned for the 03If* primary component. The bolo- 
metric albedos where set to A=1.Q (Rucinski 1969), and 
gravity brightening coefficients of value g=1.0 (Lucy 1976) 
were used. These quantities are the usual guesses for radia- 
tive envelopes. Linear limb-darkening coefficients were ob- 
tained from van Hamme (1993). A circular orbit was as- 
sumed, according to that arise from the examination of the 
light curves and is expected for such a close binary. All these 
parameters, together with the ephemerides supplied by Al- 
cock et al. (1997) (P = 1.404 740, Eq = 2 449 073.7109), 
were kept constant. 

First, the semimajor axis a, the systemic radial veloc- 
ity V-y and the mass-ratio q were adjusted using only the 
radial velocity measurements. Thereupon, these parameters 
were fixed and the orbital inclination i, the temperature of 
the secondary star T2, the luminosity of the primary Li and 
both potentials Qi and ^2 were fitted using the light curves. 
Then these two steps were repeated iteratively until the so- 
lution converged. Given that after a few iterations with the 
differential corrections program (DC), the solution evolved 
to an overcontact configuration, the program operation was 
set to mode 1. In this mode, the surface potentials of both 
stars are kept identical (Hi — ^2), the same as the gravity 
brightening coefficients {gi = 32), the bolometric albedos 
{Al — A2) and the limb-darkening coefficients {xi = 3:2). 
Besides, the polar temperature of the secondary is set by 
the gravity brightening law of the entire common envelope. 
Consequently, in the successive iterations only the parame- 
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Figure 3. Top: Observed and modelled V light curve for 
MACHO*05:34:41. 3 -69:31:39 . Bottom: (O-C) residuals for the 
light curve. 



Table 2. Model Parameters 



Ti 


50000 K (adopted) 


Si: 92 


1.00 (adopted) 


Ai,A2 


1.00 (adopted) 


a;i,a;2 


0.616 (adopted) 


e 


0.00 (adopted) 


a 


22.2 ± 0.22 Rq 


y-i 


263.1 ±3kms-l 


q (Ma /Ml) 


0.64 ±0.01 


i 


67 ±1° 


01,^2 


3.046 ± 0.03 



T2 



49490 K 



ters i, Li and Q.i (together with a, V~f and q) were fitted. 
The iterations were stopped when the corrections supplied 
by the DC program were smaller than their own errors. Fig. 
2 and Fig. 3 show the computed V light and radial velocity 
curves, together with the observational data. 

The adopted and fitted model parameters and the re- 
sulting stellar dimensions are listed in Tables 2 and 3, respec- 
tively. Fig. 4 displays the aspect of the system at different 
phases. 




2.3 Error Estimates 

The estimation of the parameters' uncertainties was at- 
tempted in two ways: First, a last run allowing the DC 
program to fit all the free parameters (i.e., a, V-,, i, fli, 
q and Li) at one time was performed, and then the stan- 
dard deviations of the diflferential corrections supplied by 

Table 3. Star Dimensions 



Ml 


41.16 ±1.2 Mq 


Ri 


9.56 ± 0.02 Rq 


Mi,o; 1 


-9.49 


loggi[c3s] 


4.09 ± 0.02 


M2 


27.04 ±1.2 Mq 


R2 


7.99 ± 0.05Rq 


Mtol 2 


-9.05 


log92[c3s] 


4.07 ±0.02 



Figure 4. The aspect of the MACHO*05:34:41.3-69:31:39 over- 
contact system. Top: at primary minimum, middle: at first 
quadrature and bottom: at secondary minimum. 



the DC program were used as error estimates. Second, the 
differences between two independent fits, performed using 
the radial velocity data together with only the V or R light 
curve separately, was used as error estimate. These differ- 
ences were smaller than the errors estimated by the first 
method, which are quoted in tables 2 and 3. 

Weighing the uncertainties of the photometric and spec- 
troscopic data, these error estimates seem rather low. This is 
not surprising, since in operation mode 1 many parameter 
constraints are applied, leaving only the six above quoted 
quantities to be adjusted. Taking into account this fact, our 
resulting parameters (and mainly our error estimates) must 
be considered with caution, since this very close system, with 
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a ~ 1.4 day orbital period and an 03 star, surely presents 
appreciable departures from hydrostatic equilibrium geom- 
etry due to the strong stellar winds and radiation pressure. 

Since this star probably belongs to an obscured tight 
star cluster (the resolution of the Digitized Sky Survey and 
the MACHO plates does not allow to address this matter), 
also it was attempted to include a third light in modelling 
the hght curves. For the best fit, the third light represents 
only ~ 2.6% of the fiux at quadratures, and the other pa- 
rameters of the model do not change meaningly. 

Looking at the (O-C) residuals in Fig. 2, two slight sys- 
tematic trends can be observed. On one hand, the secondary 
minimum seems not as deep as it is modelled by the Wilson- 
Devinney program. This feature indicates a bigger tempera- 
ture difference between the components, in accord with the 
spectral types. The mode of operation 3 of the DC pro- 
gram retains all the parameter constraints for overcontact 
systems, except for the secondary's temperature T2. An at- 
tempt to estimate the secondary temperature from the light 
curves yielded T2 = 47200 ± 1700 K (for Ti = 50000 K). 
However, since the difference in surface temperatures affects 
mainly the far UV brightness in such very hot stars (in which 
it is urgent to acquire data) and its influence in the V and R 
bands is slight, this temperature estimation may be consid- 
ered with caution. Taken into account the photometric er- 
rors, a temperature difference greater than 5000 K between 
the two components can be rejected. The other systematic 
deviation of our modelled light curve (very slight also) is an 
excess of brightness in the maximum following the secondary 
minimum (O'Connell effect, see Davidge & Milone 1984). 
The O'Connell effect is often modelled adding "hot spots" 
on the surface of the stars, although its physical causes in 
very hot binaries are not clear. It was not attempted to be 
modelled this feature, since it is also hardly noticeable in 
comparison with the photometric errors. 



tion, especially UV photometry and high dispersion spec- 
troscopy. 
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3 RESULTS 

It was an absolutely unexpected discovery to find an 03 ir- 
type star in such a very close binary pair. There are not 
published evolutionary models suitable for this system, and 
single-star models are not adequate. 

The derived masses (~ 41 and ~ 27 Mq, with an incli- 
nation of 67°) seem to be surprisingly low. 

Alcock et al. (1997) presented a fit of the light curves 
performed by means of the Nelson-Davis-Etzel model (Nel- 
son & Davis 1972, Popper & Etzel 1981). Assuming a mass- 
ratio q = 0.95, they derived an inclination of ~ 57°. The 
difference with the results presented here is not surprising, 
given that the Nelson-Davis-Etzel model is not adequate for 
highly distorted systems and the value of q derived from the 
radial velocity measurements is significantly different from 
the one assumed in Alcock et al. 

The system resembles the 04f+06V binary Sk-67°105 
(see Niemela & Morrell 1986, Haefner et al. 1994) but it is 
hotter and closer. On the other hand, our mass estimate for 
the primary is similar to that derived by Antokhina et al. 
(2000) for the 03 star of the non-eclipsing system HD 93205 
(~ 45 Mq for their best fit). 

Undoubtedly, this system will require further investiga- 
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